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We present physical arguments as well as numerical evidence to show that the Cooper pair l^gBp) 
function vanishes on the Fermi surface in strongly correlated electron systems such as t — J and ka^e- 
U Hubbard models in one and two dimensions. Using exact diagonalization, we study a correlation 
function which is a measure of the gap function even in a finite system, and present numerical results 
for this correlation function on small clusters of 8 (Hubbard) and 16 {t — J) sites. 



The study of high-Tc cuprates has brought to light a class of singlet superconductors with properties rather different 
from the BCS type of superconductor. In this context, there have been intensive photoemission, muon spin resonance 
and tunneling studies of high- Tc superconductors in recent times, as well as several theoretical suggestions. [1] In 
particular, some of the recent models of high- Tc superconductors have given rise to the possibility of gap functions 

Ak = {cl.cl^,) (1) 

having unusual symmetry properties, e.g. Ak which is odd across the Fermi surface (FS), [2,4] although the issue of 
the symmetry of the gap function still remains unresolved. In the context of the Resonating Valence Bond (RVB) 
theory, which is based on strong electron correlations, one has always looked for unique signatures in some of the 
physical properties. The behaviour of the superconducting gap function in k-space offers one such signature. 

Anderson [4] has pointed out in several contexts that the gap function vanishes on the FS and has the property 
sgnACkp + (5k) = —sgnACkp — (5k). At simple mean-field level, this is true for the t ~ J model only at half filling 
(Mott insulator). For the doped case, the analysis is difficult and Anderson did not provide a convincing argument 
as to why the gap should vanish on the FS away from half-filling. [7] We must point out here that the possibility of 
odd behavior as a function of (|k| — kj?) was anticipated by Cohen [3] many years ago, and recently many authors [2] 
have looked at a more general type of singlet pairing that is odd in k and to. 

The aim of the present letter is to provide physical arguments and numerical evidence for this claim. Wc point 
out that the singular forward scattering between two electrons with opposite spins on the FS, arising from strong 
correlation, has an interesting consequence on the structure of the gap function around the FS; in particular, it has 
a tendency to vanish on the FS, independent of the rotational symmetry in k-space. Our numerical studies on small 
clusters yield results which are indeed consistent with these expectations. [8] These results have been obtained through 
exact diagonalization on small clusters in ID (8-site Hubbard and 16-site t — J chains) and in 2D {^/S x -s/S Hubbard 
and 4 X 4 i — J planes) . Within the finite-size limitation, we clearly see a change of sign of the gap functioifliailong 
radial directions across the FS, indicating the presence of a line of zeros close to the FS, and furthermore, it has 
s-wave symmetry. The key argument of this paper is that the nature of charge fluctuations in k-space imposed by real 
space repulsions is incompatible with the nature of charge fluctuations demanded by a finite Ak on the FS. 

Our starting point is the one-band large-C/ Hubbard model or the t — J model, which are believed to describe the 
low-energy physics of high- Tc superconductors. Let us begin by noting that in the large-?/ Hubbard model, the on-site 
pairing amplitude A^^ = i'^i '^li) diminished considerably on account of the restriction on double occupancy. In the 
t — J model, of course, double o. 



A., = {cl^cl)=G ■ (2) 



In general, there is no reason to expect Ak to vanish on the Fermi surface. In a simple d-wave superconductor in 
a 2D square lattice, for example, the gap goes to zero along the lines = ky and k^ = —ky and these lines intersect 
the FS at four points. But, is there any special advantage for the line of zeros to coincide with the FS in 2D and in 
ID? We will now argue that the answer is in the affirmative. 

Let us begin by examining the behaviour of the BCS wave function close to the FS: 

\BCS)=l[iu^ + Vkclcl^^)\0) (4) 

k 

which may be rewritten as 

\BCS) = n '{< + ^/2TO&Lk + vlbl-A.-i.m, 

k 

where the product is now over only half of the k-space (e.g. fc^^ > and all ky in 2D), with 

^L-k = ^(clcTC-ki - 4icLkT) (6) 

as the singlet pair creation operator on points k and — k in k-space. It is clear that is the proMlHility amplitude of 
finding no singlet pair with momentum (k, — k), u^vu is that of finding one singlet pair (of charge 2e), and that of 
finding two singlet pairs (of total charge 4e). The BCS state has identical phase relations for various configurations 
of pair occupancy in k-space. That is, when the product in Eq. (5) is expanded out, the resulting sum has identical 
phase for all terms, each term corresponding to different configurations of the (k, — k) occupancy. Superconductivity 
can thus be thought of as a coherent charge-2e fluctuating state in k-space. Since Uk^^k is non-zero only in a thin 
energy shell around the FS, the coherent 2e charge fluctuation is concentrated around the FS. (It is interesting to 
note that this coherence in k-space results in phase coherence among the Cooper pairs in real space also). Away from 
the shell, we either have a completely filled band (inside the FS) or a completely empty band (outside the FS) and 
hence no charge fluctuations. Thus, Ukfk(~ Ak) is a measure of coherent charge fluctuations in k-space. 

We now argue that strong correlations in real space lead to a suppression of such coherent charge fluctuations close 
to the FS. Strongly correlated electrons in ID and 2D, described by a large-JJ Hubbard model, have certain unique 
features close to the FS. It is well known that in the ID Hubbard model there is singular forward scattering between 
two electrons with opposite spins close to the FS. This leads to a finite phase shift [5] at the FS and the consequent 
failure of the Fermi liquid theory, resulting in the vanishing of the discontinuity in nk at the FS (Luttinger liquid 
behaviour). It also implies an effective hard-core repulsive pseudopotential between electrons with opposite spins 
close to the FS. Thus no two electrons close to the FS, with opposite spins, can have the same momentum, thereby 
making k-points close to the FS essentially singly occupied. Coherent pair fluctuations are thus unlikely to develop 
on or very close to the FS, but are not forbidden away from the FS. It is therefore likely that the line of zeros of Ak 
implied by the global constraint (Eq. (3)) will coincide with the FS. 

Generalizations of the above argument to two and higher dimensions is straightforward provided the following is 
true: the projective consraint of no double occupancy in real space should lead to singular forward scattering and 
consequent failure of Fermi liquid theory. As discussed earlier, the single occupancy constraint around the FS in 
k-space is likely to force the (d — 1) dimension nodal surface of Ak to coincide with the {d — 1) dimensional FS. In 
particular, the 2-d case can be understood in the spirit of Anderson's tomographic Luttinger liquid picture [6] where 
we have a collection of 1-d chains in k-space. 

The above discussion can be quantified by studying the behaviour of Ak numerically, for which purpose we perform 
exact diagonalization calculations on small clusters of 8 and 16 sites. However, it is not possible for us to calculate 
Ak directly, since we have chosen to work in a number-conserving basis. We thus define a correlation function Ak,k' 
as 



(5) 



Thus if Ak vanishes on the FS and is an odd function across the FS, it can be seen that Ak,k' > when k, k' both 
he either completely inside or completely outside the FS, Ak,k' = when either of k, k' lie on the FS, and Ak,k' < 
when one is inside and the other is outside. In a strongly correlated system, we do not expect a complete factorization 
as in Eq. (8), but a form 

Ak,k' - AkA;;, +gk.k', (9) 

where (7k, k' is expected to be small when k and k' are well separated. We shall therefore infer the behaviour of Ak 
from exact results on Ak,k' by assuming a factorization of the form above. 

We have evaluated Ak,k' using exact diagonalization on 8-site Hubbard and 16-site t — J chains in ID and -s/S x \/8 
Hubbard and A x A t — J planes in 2D, for all k, k' within the Brillouin zone. Assuming that Ak changes sign across 
the FS and a factorization of the form Eq.(8) we expect to see a change in the sign of Ak,k' as k is varied across the 
FS for each value of k'. We present results for a fixed value of k' = and vary k along appropriate direction(s) across 
the FS. The k' 7^ results are qualitatively similar. 

The results for ID are particularly striking. Figure 1 depicts Ak,o for an 8-site Hubbard chain with U — 10 for 0, 
2, 4 and 6 holes. A clear sign change is observed across kp. To contrast this with the U < Hubbard model, we 
have also plotted Ak.o for U = —10 with 2 holes in the same figure. As expected, Ak,k' is always positive and nearly 
flat around k^?. This is consistent with on-site pairing for ^ >> 1, which implies (c|^|cj|) ~ Sij making Ak virtually 
k-independent. Figure 2 depicts Ak,o for a 16-site t — J chain with 0, 2 and 12 holes for J — 0.32, which again shows 
a clear change of sign across k^- Figure 3 is a plot of kg vs. kp for the 1-d chains, where kg is the node of Ak. An 
approximate linear relationship is clearly evident, with the spread being due to the degeneracy of the ground state. 
In two dimensions, we present results for the 4x4 cluster in Fig. 4 and 5. 

Figure 4 shows Ak,o for the 4 x 4 t - J cluster (J = 0.32), again along (0, 0) - (f , f ) - (tt, tt), for 0, 2 and 12 holes. 
As before, Ak,o changes sign across hp for all the three cases. We have observed a similar sign change along the 
direction (0, 0) — (|-, 0) — (tt, 0). Figure 5 is a surface plot of Ak,o for all k = {k^, ky), — tt < k^, ky < tt. In addition 
to the sign change, this also indicates the angular symmetry of Ak . As k is moved along a closed path around the 
origin, we do not see any sign change in Ak,o , which clearly indicates that the order parameter has s-wave symmetry. 
One consequence of this is that d-symmetry need no longer be invokeH-to enforce A^ = 0. The system prefers to 
have s-symmetry and the sign change across the FS satifies the global constraint. 

To summarise, we have investigated the behaviour of the gap function Ak across the Fermi surface in the t — J 
and large-?/ Hubbard models in one and two dimensions, at and away from half-filling. We have presented physical 
arguments for the vanishing of Ak on the FS. Exact diagonalization results on 8- and 16-site clusters, both in ID and 
2D, clearly indicate a change of sign of Ak across kj?, and in addition, the 2D results have s-wave symmetry. Our 
results, which essentially gives a gapless s-wave superconductor, has consequences for several experimental results in 
the cuprate superconductors below - in particular, the temperature dependence of the London penetration depth, 
NMR relaxation rate, zero bias conductivity in tunneling experiments, specific heat etc. Some of these consequences 
are being investigated. 
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Figure Captions 

1. Ak,o as a function of ^ for an 8-site Hubbard chain {U = 10) for 0, 2, 4 and 6 iioles. The corresponding 
kp = and respectively. For comparison, we also plot Ak,o for C/ = —10 for 2 holes. 

2. Ak^o as a function of ^ for an 16-site t— J chain (J = 0.32) for 0, 2 and 12 holes. The corresponding = ^, ^ 
and I respectively. 

3. The kg vs. kp curve for the 1-d chains; ks is the node of Ak- The straight line is the kg = kp plot as a "guide 
to the eye" . 

4. Ak,o as a function of ^ for a 2D 4 x 4 t — J cluster (J = 0.32) for 0, 2 and 12 holes, k is varied along the 
direction {x,x), with a; = 0, |,7r. \s.p corresponds to x = ^ for and 2 holes. For 12 holes, (0,0) is completely 
inside the FS and the other two points completely outside. 

5. A surface plot of Ak,o as a function of ^ for the 4 x 4 t — J cluster (J — 0.32) for 2 holes, kp corresponds 
to the set of points (±f , ±f ), (±|, 0), (0, ±f ). 



